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1. Introduction 

In the 1970's and early 1980's, there was a significant amount of interest in the 
radiative lifetimes of Rydberg states of alkali metal atoms. From the experimental 
perspective, alkali atoms are relatively easy to work with. High atomic densities can be 
achieved easily in a cell or atomic beam, and the wavelengths needed to photo-excite 
specific Rydberg states coincide well with those available from the (then) newly-invented 
tunable, pulsed dye lasers [H; 0]. The state of interest (termed the "target state" in this 
paper) could thus be excited with a temporal precision of 10 ns or less, and the radiative 
decay could be monitored by observing the emitted fluorescence using a photomultiplier. 
There is theoretical interest in radiative lifetime measurements of Rydberg states of 
alkalis since their wave functions are very similar to those of hydrogen, and calculation 
of atomic properties of alkalis is much less computationally intensive than for other 
groups of elements. Experimental values for radiative lifetimes provide useful data for 
testing the accuracy of the wave functions used in the calculations 

A comprehensive list of Rydberg state lifetime measurements in Li, Na, K, Rb, 
and Cs as it existed in 1984 can be found in Theodosiou's theoretical paper [j]. A 
very significant issue that became apparent as higher n-states (where n is the principal 
quantum number) were measured was the effect of black body radiation (BBR) on 
Rydberg state lifetimes. BBR has a negligible effect on low-lying energy states, whose 
radiative lifetimes are dominated by spontaneous electric dipole transitions to the lowest 
atomic states allowed by the selection rules. However, dipole matrix elements for n,i ^ 
n±l,i±l transitions increase approximately as n^. In addition, the energy spacing of 
states decreases as n~^, and thus the photon energy needed to drive transitions between 
Rydberg states tends towards the peak of the BBR spectrum as n increases. Above 
n ^ 10, therefore, BBR stimulated depopulation becomes non-negligible compared 
to the spontaneous decay of a state. Experimentally, the effect of BBR became 
apparent when the measured Rydberg state lifetimes were found to be significantly 
lower than those calculated assuming purely spontaneous decay, i.e., ignoring stimulated 
depopulation by BBR. The ability of BBR to significantly affect Rydberg state lifetimes 
of alkalis was investigated in a number of experimental studies in the late 1970's 0; S l^l. 
In particular. Spencer et al made a thorough study of the effect of changing the BBR 



spectral temperature on the state lifetime 10|. Calculations of BBR redistribution 



rates, and the effect of BBR on Rydberg state lifetimes, were made by Farley and Wing 



11 



and Theodosiou 0]. Additionally, Galvez et al fl^ compared the experimentally 
observed population redistribution due to BBR stimulated transitions with calculated 
distributions based on theoretical rate values obtained by similar methods to those 
reported in 0] and 11 . 

Measurements of Rydberg state lifetimes that use thermal samples are limited to 
states below n ~ 20. This is because the radiative lifetimes increase with n, and at 
n = 20, they are of order 10 fis. During this time, thermal atoms travel distances on the 
order of 1 cm. Accurate measurements require observation for several lifetimes, which 
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necessitates a significant detection volume. Additionally, collisional depopulation effects 
and superradiance become significant at long radiative lifetimes. 

In the past 10-15 years, there has been great interest in cold Rydberg atoms. When 
made in sufficiently dense samples, they exhibit many-body collision behavior 13|; Il4l |. 
as well as suppressed laser excitation probability due to an interaction-induced blockade 
(for instance the van der Waals blockade, or the dipole blockade, or some other density- 
dependent excitation suppression mechanism) 15|; [ly; [iTI; ll8|] . In addition, novel long- 
range molecules can be created 120], and dense, cold Rydberg samples have also been 
found to evolve spontaneously into ultra-cold, strongly coupled plasma 2l|; |22|; |23|. All 
of these studies start with cold atoms in a magneto-optical trap (MOT), and such a 
source is ideal for measurements of Ryd berg atorn properties such as quantum defects 



;orn_pr 

I m 



24; 25 and radiative lifetimes 26; 27; 28 



3ll; l3^. In the ultra-cold MOT 
environment, the atomic velocities are less than 1 m/s and hence the atoms remain in 
the detection volume long enough for an accurate determination of the properties of 
Rydberg states to be made. 

While various atoms have been trapped in a MOT successfully, Rydberg lifetime 



measurements using cold atoms have so far been reported only for Rb 26 



27; 28; 29; 30; 



31 and Cs 



32| . In all of these experiments except 30||, the technique that was used to 



monitor the decay of the Rydberg state population was selective field ionization (SFI) 
with a ~ 1 /xs rise-time high- voltage pulse. In 301], the technique for monitoring the 
Rydberg population was to detect photo-electrons or photo-ions created when a mid- 
infrared (~ 10 /zm wavelength) laser pulse illuminated the sample at a well-determined 
time after the target state was excited. (These techniques are discussed more fully below 
in section [3m ) These recent experiments have stimulated renewed interest in extending 
theoretical calculations to the previously unstudied region n > 30 [6]. Additionally, 
lifetime measurements have been used to distinguish exotic Rydberg molecular states 



from their parent atomic states 20 



The SFI detection technique that was used to monitor the Rydberg state population 



in [26|; l27|; l28|; l29|; l3ll; 122] must be applied carefully to avoid significant systematic errors 
in the lifetime measurement. (Photoionization by an infrared laser pulse, the technique 
used in |30|], may also lead to systematic errors.) Specifically, the arrival times of 
electrons (or ions) created by SFI at the detector must be carefully gated so that only 
those from field-ionization of the target state are observed. Due to BBR, states with 
similar energies to the target are populated by stimulated transitions during the decay 
time of the target, and if the SFI signal from these nearby states cannot be rejected, 
the method yields an inaccurate lifetime measurement. One of us has commented on 
this issue |33| with regard to the results of a group at Universidade de Sao Paulo, which 
are presented in 
implications 



26; 27; 28 



and the comment describes fully the problem and its 
In their response |34|. the Sao Paulo group clarifies their experimental 



34| . systematic uncertainty when the SFI 



procedure, but as they observe in [29( and 
detection technique is used is still possible. 

The present paper presents Rydberg state lifetime measurements of Rb for most 
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of the same states previously reported in |26|; |27|; l28|; l29| . However, our measurements 



were made with a modification of the SFI technique that allowed the decay of the 
target state to be unambiguously distinguished from signals arising from the effects of 
BBR stimulated transfer of population to other Rydberg states which also contribute to 
the SFI signal. This was achieved by using frequency-tuned pulses of millimetre-waves 
(mm-waves) to sample the population of the target state as a function of delay after the 



24 



25 has determined 



exciting laser. Previous work by a group at University of Virginia 
the frequencies of ni — > n'i' mm-wave transitions in Rb with great precision, and we 
use this knowledge to accurately identify which part of the SFI signal is characteristic 
of the target state population. 

Below, we describe our experiment and our results. First, however, we discuss 
the present state of lifetime calculations. Theoretical considerations have an important 
impact on this experiment, since our lifetime measurements are made in a BBR field 
with a spectral temperature of nominally 300 K. While our experiment (as with the 



previous cold-atom lifetime studies [26|; l27|; l28|; l29|; l30|; ImI; l32|) has no ability to 



parse the zero-Kelvin lifetime from BBR-stimulated depopulation rates, the theoretical 
calculations do distinguish these effects. Comparison of experiment with theory then 
becomes dependent on theoretical calculations for the BBR depopulation rates, and as 
we describe, the works described in 26|; |27|; |28|; |29|; all use an inaccurate equation 
for these rates. 



2. Review of previous theoretical results for Rb Rydberg state lifetimes 

The total radiative decay rate of a Rydberg state at temperature T, rgg-(T), is related to 
the radiative lifetime of that state, r^g-, by the relation rgg-(T) = 1 / t^q. It is convenient 
to consider two radiative depopulation effects that constitute the total radiative decay 
rate separately. The first effect is the decay rate in the absence of black body radiation, 
Fq = I/tq, where tq will be referred to throughout this paper as the "zero-Kelvin 
lifetime". The second effect is the stimulated black body depopulation rate due to 
radiation with a characteristic temperature T, rBB]^(T). The total radiative decay rate 
is the sum of these two effects: 

reff(T) = ro + rBBR(T). (i) 

The zero-Kelvin vacuum radiation field that produces Fq generally results in transitions 
from the target state to the lowest energy state accessible by an electric dipole transition. 
On the other hand, black body radiation stimulates electric dipole transitions from the 
target to both higher and lower energy states, and the transitions with the highest 
probability are those with the smallest change in n. 

2.1. Calculations of zero-Kelvin lifetimes 

There have been a number of Rydberg state lifetime calculations for alkali metals, but 
until very recently, there was no theoretical analysis that was directly relevant to the 
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present work on Rb. There are three principal reasons for this. First, the theoretical 
studies were confined to \ow-n values, for instance Theodosiou's work was limited 
to n < 20 for Rb, and the study of Gounand (3| considered only n < 28, though that 
of He et al. [5|] extended up to n = 50. (All three of these theory papers calculated 
lifetimes for all Group I elements except H and Fr.) Second, while Theodosiou calculated 
lifetimes for K, 350 K, and 460 K black body radiation fields, both Gounand and He 
et al calculated only tq, i.e., they ignored stimulated depopulation by BBR. Finally, 
none of these three theoretical works discussed how their results could be extended to 
arbitrary-temperature BBR fields. Very recently, however, a group in Novosibirsk [gI has 
calculated both the zero-Kelvin lifetimes and black body depopulation rates (discussed 
below) in all the alkalis at arbitrary temperature for 10 < n < 80 (their paper gives 
an extensive discussion of the calculation techniques used in previous studies which we 
omit here). 

It is usual to describe theoretical zero-Kelvin lifetimes using a simple analytical 
expression of the form 

To = Tsn^ff, (2) 

where and e are constants found by fitting the calculated Tq values as a function of 
the effective principal quantum number, rigg- = n — 6e, where Se is the quantum defect 
of the state in question [35]. Values for Tg and e reported in jl], and [6| are given in 
table [H However, fitting different ranges of UqQ yields slightly different Xg and e values. 
This led He et al to specify the maximum error for the use of the constants outside the 
range of validity given in table [1] (10 < n < 30), and also to propose an alternative 
fitting function for the behavior of tq versus rigg- 

To = ao + cti^eff + '^2^eff + (3) 

where Oq, ai, 02, and 03 are obtained from the fit and are given in table 9 of [sf. The 
maximum error when using outside the specified range of 10 < n < 30 is significantly 
lower than for ([2]). 



2.2. Calculations of black body depopulation rates 



The BBR depopulation rate, rBBj^(T) has been calculated at different levels of precision. 



Cooke and Gallagher 37| and others 26|; |27 
approximation: 



28l : I29I : |32| | have used the following 



-pApp 
^ BBR 



(4) 



where atomic units (A.U.) are assumed, and the symbols have their standard meanings. 
While the accuracy of the approximation improves as ngg- increases, the work of Farley 
and Wing ll|] reveals that it overestimates the black body decay rates by ~20% for ns 
states, ~30% for np states, and ~40% for nd states in Rb at n = 30. Unfortunately, 
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Table 1. Values of the parameters ts and e in 12]) obtained from the results of three 
major theoretical lifetime calculations of tq for Rb Rydberg states. Note that Q and 
Q do not consider the different j-values of the np and nd states, while Q calculated 
the lifetimes and fit parameters for nsi/2, npi/2, np3/2, nd^/2 and states. For the 

np3/2 states, two sets of values for Tg and e are given for the work of Beterov a/ 0]. 
The values identified by an asterisk (*) are as published in Q, while those identified 
by a dagger (|) are corrected values, and are to appear as an erratum [sBl (see section 

11231). 



Reference 


State 


Ts 

(ns) 


e 


Validity range 


Gounand [3] 


nsi/2 


1.43 


2.94 


10 < n < 28 




np 


2.76 


3.02 






nd 


2.09 


2.85 




He et al [5j 


nsi/2 


1.25 


2.99 


10 < n < 30 




np 


2.65 


2.95 






nd 


1.65 


2.84 




Beterov et al [6] 


nsi/2 


1.368 


3.0008 


10 < n < 80 




npi/2 


2.4360 


2.9989 






np3/2 


2.5341* 


3.0019* 






np3/2 


2.2214t 


3.0026'f 






nd3/2 


1.0761 


2.9898 






nd^l2 


1.0687 


2.9897 





Farley and Wing's work on BBR decay rates of hydrogen, helium, and the alkalis, 
extends only up to n = 30. 

Recently, Beterov et al [6|] calculated rBBj^(T) numerically for alkah atom Rydberg 
states in the range 10 < < 80 at several temperatures, and found a semi-empirical 
analytical function that accurately fits their results in the form 

yBet ,rj.. ^ A 2.14 X 10^° 1 

BBrU ) exp[315780 x ^/(n^gT)] - 1 ^ ^' ^ ' 

where the constants A, 5, C, and D are found from the fit and are given in table I in 

i. 



3. Experimental 

3.1. General principle of the experiment 

The goal of the experiment is to measure the lifetimes of specific target Rydberg states. 
To do this, one must monitor the population of the target state as a function of time 
delay after excitation. Additionally, the target state excitation needs to happen on a 
time scale that is short relative its radiative lifetime. (Figure [1] illustrates the excitation 
process, and the techniques that have been used to monitor the decay of the target state 
population, which are described below.) 
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Figure 1. Diagram illustrating the various techniques that have been used to measure 
Rydberg state lifetimes. The experiment begins when a target state (TS) is excited by 
a short laser pulse (PL) at t = 0. The target state population may then be monitored 
by SFI with a voltage pulse at a known time delay after t — 0, and detecting the 
resulting ions or electrons. Alternatively fluorescence (FL) from radiative decay of the 
target state can be observed. Another method that has been used is to photo-ionize 
the target state with a short infrared laser pulse (IR), again detecting the ions or 
electrons so created. Transitions stimulated by black body photons (BBR) populate 
states close in energy to the target during the delay between excitation and detection 
(see text for the consequences of BBR on lifetime measurements). Our experiment 
uses a frequency-tuned mm-wave pulse (MW) to drive a fixed fraction of the target 
state population to a higher "monitor state" (MS), which is then field ionized. (The 
atom ionization potential is represented by the line labelled IP.) 



In most Rydberg state lifetime measurements, excitation is done using a pulsed 
laser (<10 ns pulse length) to excite the target state, and this is the method which 
we use. Three methods have been used to measure the target state population time- 



dependence. Specificall 
state radiative decays [38 
pulse 



they are: (a) detection of fluorescence from one of the target 



26; 



23; l28|;l2g; 



31 



39 



4G; 



41 



42l |; (b) selective field ionization using a voltage 



32l |; and (c) photoionization using a mid-infrared laser pulse 



3G 



All three of the detection methods described above are in some way challenged when 
high-ra Rydberg states are being measured, due to population redistribution caused by 
BBR. Depopulation of the target state by BBR contributes to the state lifetime at any 
temperature above zero Kelvin; however, BBR selectively populates states that are close 
in energy to the target [37|, and the method by which the target state population is 
being monitored needs to reject the signal from these nearby states. For fluorescence 
detection, a spectrometer may pass only photons originating from the target state. The 
necessary spectrometer resolution increases for high ra-states since the energy spacing 
between adjacent states becomes smaller and thus the photon energies from states with 
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{n — ni, {n + l)i (where ni is the target) become hard to distinguish. In principle, 
SFI with a ~ 1 /is rise-time voltage pulse may be used to distinguish nearby states, 
provided that electrons (and not ions) are detected. In atomic units, the field required 



to ionize a state with effective principal quantum number n^Q is l/{lQn'^g) [35|. The 
time-of-fiight of the electrons to the detector is much less than 1 /xs, and so they arrive 
at the detector at essentially the same time that the electric field of the voltage pulse 
becomes large enough to adiabatically ionize the Rydberg state. Electrons from field 
ionization of states with different n^Q arrive at the detector at different times, and the 
original state can, in principle, be determined. However, states with similar n^Q but 



different i have different field ionization threshold behaviours 43|, and unambiguous 



correlation of electrons arriving at the detector at a certain time with a specific ni state 



is not possible at high n SSj. Finally, an intense mid-infrared laser pulse with frequency 



u (e.g., from the 10.6 fim CO2 laser used in [30t|) will photoionize almost any state closer 



to the ionization limit than hu. This technique therefore has no state specificity. 

The detection strategy of our experiment is shown schematically in figure [1], and 
is a modification of the SFI detection technique. We avoid the limitation in the SFI 
method by applying a short (~ 1 fis) pulse of mm-waves just before the SFI pulse 
to "tag" the target state. We use this ability to identify the part of the SFI electron 
signal that unambiguously mirrors the target state population, and rejects signal from 
nearby states populated by BBR. Briefiy, a short (0.5 - 5 fis), frequency tuned, mm- 
wave pulse is used to excite a fixed fraction of the target state population to a higher 
"monitor" state at a well-determined delay after the target state is populated. We then 
field ionize the monitor state with a voltage pulse at a fixed time after the end of the 
mm-wave pulse. The SFI pulse amplitude is set so that the monitor state is ionized, but 
not the target, and we gate and detect the electron signal at a specific time-of-fiight. 
As described below, these two techniques (mm-wave excitation and gating the electron 
signal) provide the necessary discrimination against signals originating from states close 
in energy to the target. 

We now describe the four main parts of the apparatus. Specifically, these are the 
magneto-optical trap, the pulsed laser system, the mm-wave apparatus, and the timing 
and detection electronics. 

3.2. Magneto- optical trap 

Our apparatus is basically a vapour-cell MOT [4^] with internal field plates for applying 



electric fields and a micro-channel plate detector (MCP) [iSj; . A vapour of rubidium 
atoms is contained in a cylindrical stainless steel vacuum chamber of diameter 150 mm 
and height 200 mm. The chamber is pumped by a 20 1/s ion pump to a base pressure of 
less than 1 x 10^^ torr. There is ample optical access for the trapping laser beams, as 
well as for pulsed lasers and mm-wave radiation. There is a magnetic field gradient of up 
to 0.2 T/m at the centre of the chamber, created by a pair of air-cooled anti-Helmholtz 
coils outside the chamber. The source of rubidium atoms is a side-arm heated to 60 - 



Radiative lifetime measurements of rubidium Rydberg states 



9 



80° C. 

The Rb atoms are cooled by three orthogonal, retro-reflected laser beams of the 
appropriate circular polarizations. These beams are derived from a single external cavity 
diode laser [45| that is locked to the appropriate hyperfine component of the 5si/2 — >■ 
5p3/2 transition in ^^Rb. This light is amplified (after suitable isolation) by an injection- 
locked non-cavity diode laser, expanded to a diameter of 10 mm, and split into the 
three orthogonal beams that are introduced into the vacuum chamber. Each trapping 
laser beam has a power of ~ 20 mW. A second amplified external-cavity laser system 
is used as a repumper. The trapped atoms are observed using a video camera, and we 
measure the diameter of the cloud of atoms by imaging the fluorescence from the atoms 
onto a linear diode array. We determine the number of trapped atoms by measuring 
the fluorescence emitted into a known solid angle using a sensitive optical power meter. 
Using this apparatus, we trap up to 1.2 x 10^ ^^Rb atoms in a spherical volume of 
diameter fal.5 mm (FWHM inferred from diode array signal). The maximum cold-atom 
number density is therefore ~ 5 x 10^° cm~'^, though we typically run our experiments 
at number densities which are at least ten times lower than this to avoid collisional 
depopulation of the target state. Rubidium atoms trapped by this method typically 
have temperatures of ~ 100 /iK, corresponding to atom velocities of < 0.2 m/s. 

The atoms are trapped midway between two parallel field plates made from high- 
transparency stainless steel mesh, separated by 18.3 mm. Excitation of 5p3/2 atoms to 
a Rydberg state is achieved using a tunable pulsed laser (see section [3731) . and this state 
is field ionized at a specified delay after excitation using a ~ 1000 V, 1 /xs rise time 
voltage pulse applied to the parallel meshes. The resulting electrons are detected using 
the MCP. 



3.3. Pulsed laser 

We use a blue (^480 nm) pulsed laser to excite ns and nd Rydberg states of '^^Rb starting 
from the 5p ^^3/2 state. The np states are also accessible provided a small (1-5 V/cm) 
1 /is-long electric field pulse is applied simultaneously with the laser pulse. (The 5p3/2 
state is the upper state of the MOT trapping and repump transitions, and close to 50% 
of the cold Rb atoms in the MOT are in this state.) For 5p3/2 — > ni transitions, the 
Doppler widths of optical transitions for 100 /iK Rb atoms are less than 1 MHz. The 
spectral width of such a transition is thus limited by the 5p3/2 state's 5 MHz natural 
width, since the ni Rydberg states have much longer lifetimes than the 5ps/2 state. For 
this reason, many groups excite cold Rydberg states using continuous- wave (cw) lasers 



16 



17 



18 



20 



because the 1 MHz linewidth available from commercial cw dye, diode, 
and titanium-sapphire lasers is ideal for exciting and performing spectroscopy on cold 
Rydberg atoms. 

On the other hand, it is difficult to make the short pulse of light (< 100 ns long) 
needed to excite the target state for a radiative lifetime measurement from a cw laser 
source. Instead, we used a variation of a Nd:YAG pumped pulsed dye laser. However, 
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Figure 2. Schematic diagram of the laser system: ECDL, 960 nm external cavity 
diode laser; 01, optical isolator; SMF, single-mode fibre; L, lens; WM, wavemeter; CE, 
confocal etalon; DC, dye cuvette with LDS 925 dye; M, mirror; PB, Pellin-Broca prism; 
G, 532 nm pump light from Nd:YAG laser; X, KNbOa frequency doubling crystal; BS, 
beam splitter. 



when exciting cold Rydberg atoms using a ~10 ns pulsed dye laser (for instance a 
Hansch- [l| or a Littmann-type j3] dye laser) shot-to-shot fluctuations in the target 
state population are very large (±100% is not unusual). This is due mainly to the 
rapidly changing spectral overlap between the longitudinal modes of the dye laser and 
the spectral width of the atomic transition from one laser shot to the next, an issue 
discussed fully in [2l|]. To reduce the shot-to-shot excitation efficiency variation, we seed 
a Nd:YAG pumped dye amplifier with light from an external cavity diode laser (ECDL). 
The output of this amplifier is frequency-doubled to produce narrow-bandwidth pulses 
of 480 nm hght H S • 

A schematic diagram of our laser system appears in figure [2l The ECDL (Sacher 
Lasertechnik TEC-520-0960-100) wavelength was tunable between 967 nm and 958 nm, 
which, when frequency doubled, can excite 5p3/2 ni where 25 < n < oo. Light from 
the laser is passed through two Faraday isolators (>72 dB of isolation), and a small 
part is split off and directed towards a wavemeter (Burleigh WA-1500) and a 0.5 meter 
confocal etalon. For optimal stability, the ECDL is located on an optical table, but 
the dye amplifier is much less sensitive to acoustic and mechanical noise caused by the 
Nd:YAG pump laser. We therefore convey the 960 nm light from the ECDL through a 
fibre to the table where the dye amplifier is located. We use three fiowing dye amplifier 
cells containing LDS925 dye. The first cell has a concentration of 0.06 g/1 of LDS 925 
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dissolved in 15% propylene carbonate and 85% ethylene glycol. The latter two cells 
have a dye concentration of 0.12 g/1. The cells are pumped by 532 nm light from a 20 
Hz Nd:YAG laser (Continuum Surelite-II), and the energy/pulse directed to each cell 
(from the first stage to the third stage) are 20 mJ, 10 mJ, and 20 mJ respectively. We 
then frequency double the amplified 960 nm light in a 3 x 3 x 5 mm^ KNbOs crystal, 
and separate out the 480 nm light using a Pellin-Broca prism. Using this apparatus and 
seeding the amplifier with 10 - 15 mW of 960 nm light (after the fibre and isolators), 
we obtain as much as 200 /iJ/pulse with new dye, but more typically 50 - 100 /iJ/pulse. 
The spectral width of the 480 nm light is Fourier transform limited by the 10 ns Nd: YAG 
pulse to approximately 200 MHz. (We can easily resolve the Rb 32^3/2 and 32^5/2 states 
which are separated by 370 MHz.) The ^480 nm light from the pulsed laser is directed 
into the centre of the MOT chamber via a 1-metre focal length lens which is defocussed 
from the centre of the trap so that all the cold atoms are illuminated. 

Using this system, the shot-to-shot variation in the Rydberg signal is < ±10%. 
(This variation could perhaps be improved by pumping the amplifier with an injection- 
seeded Nd:YAG laser, but the present stability is adequate for our purposes.) To excite 
a particular 5^3/2 ni transition, the ECDL needs to be tuned to within 100 MHz 
of the correct transition frequency. The necessar y E CDL frequencies were calculated 



using known physical data for Rb 2J; |25|; |46|; |47|; |48|; |49|| , and the laser frequency was 



measured to a precision of ±30 MHz using the wavemeter. 



3.4- The millimetre-wave system 



A schematic diagram of our mm- wave source can be seen in figure [3] [22|; |23|; l2J; |25[ . It 
can generate some 10 mW of radiation in the range 20 - 110 GHz and, at lower power 
levels, frequencies of up to 200 GHz. The source is based on a 10 MHz - 26.5 GHz 
sweep oscillator capable of producing 200 mW from 10 - 20 GHz, and 10 mW from 
20 - 26.5 GHz (Hewlett-Packard HP83630). The sweep oscillator has a resolution of 1 
Hz. The output from the sweep oscillator is switched to make pulses of arbitrary length 
(General Microwave F9114A; 0.5 - 5 /is pulse lengths are typical in this experiment), and 
then amplified and/or frequency multiplied using active or passive GaAs diode doublers 
or triplers (Pacific Millimeter and Narda) We use four different combinations 

of multipliers and amplifiers to cover the 20 - 110 GHz range, each optimized to a 
particular frequency interval. For the range 20 - 40 GHz, only an active doubler (Narda 
DBS-2640X220) is used; for 50 - 75 GHz and 75 - 110 GHz, the active doubler is 
followed by a passive doubler (Pacific Millimeter V2W0) and a passive tripler (Pacific 
Millimeter W3W0) respectively. For 40 - 54 GHz, we amplify the sweep oscillator output 
(Narda DBS-0618N223), and then use a passive tripler (Pacific Millimeter U3). Once the 
appropriate frequency is generated, the mm-waves are directed into the MOT chamber 
using a horn and gold-coated off-axis paraboloid reflector through a 150 mm diameter 
viewport. 

Various mm-wave transitions were used to populate the monitor state from the 
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MW 




PD/PT 



Figure 3. Schematic diagram of mm-wave system: MW, microwave sweep oscillator; 
CC, coaxial cable (SMA); S, microwave switch; TP, trigger pulse; I, isolator; A/AD, 
amplifier or active doubler; WG, waveguide; PD/PT, passive doubler or passive tripler; 
A, attenuator; H, horn; PR, paraboloid reflector; W, MOT vacuum chamber window; 
Rb, cold Rb atoms. 



target. For the nsi/2 lifetime measurements, we used either a nsi/2 —>■ np3/2 single- 
photon transition or a nsi/2 —>■ {n + l)si/2 two-photon transition. For the 
measurements, either np3/2 — > (n + l)si/2 (one-photon) or np3/2 — > + 1)^3/2 (two- 
photon) transitions were used. Finally, for the nd^/2 measurements, either nd^/2 
(n — 2)/7/2 (one-photon) or nd^/2 {n + l)d^i2 (two-photon) transitions were employed. 

3.5. Timing and detection electronics, data acquisition and analysis 

The 480 nm laser pulse, the mm-wave pulse, SFI pulse, and the detection electronics 
(boxcar amplifier and digital oscilloscope) are all synchronized using two digital delay 
generators (DDG; Stanford Research Systems DDG-535). One DDG is used a master 
oscillator, and is internally triggered at 20 Hz. The master DDG triggers the fiashlamp 
and Q-switch of the Nd:YAG laser used to pump the dye amplifier for the ECDL (see 
section [373]) . The master DDG also triggers an analog pulse generator which is used to 
apply a small (< 5 V/cm), 1 /xs long pulse to the Rb atoms just before the arrival of 
the 480 nm laser pulse when np states are being excited. The rest of the experiment is 
synchronized by the second DDG, which is triggered (at time t = 0) by the master DDG. 
The second DDG provides three output pulses that synchronize the mm- waves and the 
detection apparatus, as diagrammed in figure |H One pulse controls the microwave switch 
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Figure 4. Schematic diagram of timing of experiment. In (a), a partial energy diagram 
of ^^Rb is sliown, along with the excitation scheme of the target state (whose lifetime 
is to be determined): DL, MOT diode lasers (trapping and repump lasers) drive the 
5si/2 —> 5p3/2 transition; PL, pulsed laser (480 nm, 10 ns); MW, mm-wave pulse; TS, 
target state, MS, monitor state; IP, ionization potential. In (b), the timing sequence is 
shown: t, time; VP, optional voltage pulse (maximum electric field V/cm, length ssl 
/is) allowing excitation of np Rydberg states from 6^3/2! SFI, selective field ionization 
pulse ionizes monitor state but not target state; A, delay of mm-wave pulse after 480 
nm laser pulse; A, mm-wave pulse length; B delay of SFI pulse after mm- waves; all 
other symbols as in (a). 



in the mm-wave system (see section [3^ . This pulse is delayed by A after t = 0, and has 
a duration A of between 0.5 /is and 5 /is, depending on the state being measured (we 
ensure that A <^ t^q) . At a fixed delay of B = 2 fis after the microwave pulse is switched 
off, the second pulse from the DDG triggers the SFI pulse. The SFI pulse is applied 
to one of the two parallel meshes surrounding the cold atoms in the MOT (see section 
13. 2p . and ionizes the monitor state (and all higher energy states), but does not ionize 
the target state. The resulting electrons are pushed towards the MCP. The output of 
the MCP is amplified and passed to a digital oscilloscope and a boxcar amplifier. (The 
third pulse from the DDG triggers the digital oscilloscope and boxcar amplifier 1 /xs 
before the SFI pulse is triggered.) The boxcar amplifier is set to acquire the part of the 
electron time-of-flight signal corresponding to field ionization of the monitor state. 

The second DDG, the microwave sweep oscillator, and the boxcar amplifier are 
controlled by a computer via a GPIB interface using a program written in the LabVIEW 
environment (National Instruments). The program controls the microwave oscillator 
output frequency and power level, and can turn the microwave output on and off. The 
program also sets the mm-wave pulse length A, and steps the delay A from an initial 
value to a maximum value in 20 equal steps. Typically, the initial value of A is between 
1 and 5 fis, and the maximum value of A is between 60 and 200 /is later than the initial 
value, depending on the state being measured. At each delay step A, LabVIEW acquires 
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the time-gated MCP signal from the monitor state for a predetermined number of laser 
pulses (usually 10) when the microwave oscillator is on (resulting in a measurement 
S'(A)on) and when it is off (S'(A)Qff). The means and standard deviations of the means 
of 5'(A)on and S'(A)og- for the 10 laser shots are calculated. The delay is then stepped 
back to "zero" (this means the initial value of A), and the corresponding means and 
standard deviations of the means are calculated for ^(A = 0)on and S{A = O)off over 
10 laser shots each. For this A, LabVIEW calculates the fraction F 

S{A = 0)on - S{A = O)off' ^ ^ 

along with its standard deviation (ctf). Subtracting the signal acquired when the mm- 
waves are off allows us to reject spurious signals from electrons originating from SFI of 
states other than the monitor state, and from BBR stimulated population of the monitor 
state from the target state, neither of which is a constant fraction of the target state 
population. In addition, repetitively stepping A back to zero allows us to normalize 
against slow variations in the target state population, due, for instance, to changes in 
the MOT atom density and variations in the 480 nm laser pulse energy and frequency. 
LabVIEW then steps A to the next delay value, and the process is repeated until we 
have 21 F and ap values for 21 different A values. The F versus A curve is then plotted 
and fitted to find the decay lifetime using Igor Pro ( Wavemetrics) . The uncertainty in 
the A measurement is negligible due to the high precision of the DDG's used in the 
timing electronics {a a ^ 10 ns), while the uncertainty in each F value is ap, and the 
ap values are used to weight the decay curve fit. 

We measure the lifetime of each state between 5 and 40 times. From each state's set 
of lifetime measurements, we calculate the mean of the set, which is our experimental 
lifetime determination, and also the standard deviation of the mean. This standard 
deviation of the mean is combined in quadrature with the mean uncertainty arising from 
the fitting process to find the uncertainty in the lifetime measurement. The resulting 
experimental uncertainty is between 2% and 8% of the measured lifetime, with around 
5% being typical. 

4. Results and discussion 

In the following sections, we present our results, and compare them with other recent 
experimental results. We also extract the zero-Kelvin lifetimes and decay rates from 
our results using Beterov's equation for rBBj^(T) (see ([5]) above), and compare them 
with previous theoretical results. Finally, we discuss the impact of the BBR spectral 
temperature on the inferred zero-Kelvin lifetimes. 

4-1- Experimental results 

Our experimental lifetimes for the Rb nsi/2, "^^3/2 and states are listed in table El 
and plotted in figure [S] We obtained lifetimes for all the nsi/2 states in the range 
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28 < n < 45 (18 states), for states in the range 34 < n < 44 (11 states), and for 

nd5/2 states in the range 29 < n < 44 (16 states). This range of states is hmited by 
several factors. At low n, the energy spacing of adjacent levels exceeds the maximum 
mm- wave photon energy available from our apparatus (around 200 GHz), even when 
two-photon transitions are used. At hi gh n , collisions become an issue due to long 
range dipole-dipole interactions [l^; E^:l2ll: 22 : 23 : 35 1. Though we took extensive 



measures to work at low background gas pressure and low trap density, we noticed 
that the lifetime measurements tended towards a limit as n increased above 45-50, an 
effect that we attributed to collisions. (This is apparent even though our measurements 
for these states were made with the Rb reservoir valved off from the MOT chamber, 
i.e., we trapped from the residual Rb atoms in the vacuum chamber.) Another issue 
that came into play above n ~ 45 was the effect of microwaves "leaking through" the 
switch that synchronized the mm-waves with the pulsed laser. As n increases, Rydberg 
states become more and more sensitive to microwave radiation, and the small amount of 
leak-through radiation was enough to affect the measured lifetime. For nd states with 
n > 39, the bandwidth of the pulsed laser which excites the target state at t = (~ 200 
MHz) exceeds the separations of the j = 3/2 and j = 5/2 sub-levels, which is close to 
200 MHz at 72 = 39 and approximately 120 MHz at n = 44 2J] . These nd state lifetime 
measurements may therefore be subject to some averaging over the different j-levels. 
The mm-wave source that was used to populate the monitor state has a resolution that 
is more than enough to distinguish the nd3/2 —>■ {n + l)d3/2 two-photon transitions from 
the nd5/2 — (n + l)d5/2 two-photon transitions that were used for these states, which 
are separated by between 20 MHz and 10 MHz in the range 39 < n < 44. However, it 
is possible that power broadening may result in some j = 3/2 character being present 
in the nominally "pure" j = 5/2 monitor state signal for these target states. On the 
other hand, the calculations of Beterov et al indicate that the lifetimes of the two 45dj 



sub-levels differ by less than 0.2 /is at 300 K (see table VII in [6j), so the impact of any 
inadvertent averaging on our results is probably negligible. 

Two other features of our experiment have a potential to impact our lifetime 
measurements. First, our experiments were performed without switching the anti- 
Helmholtz coils off for a short time in synchronism with the laser pulse. Hence, the 
MOT magnetic field was on at all times, which could lead to state mixing due to the 
Zeeman interaction. Given the spatial size of the atoms in the MOT 1 mm), and the 
fact that the atoms were trapped at the zero-point of a magnetic field with a gradient 
of < 0.2 T/m, the maximum field magnitude experience by the Rydberg atoms was 
^ ]^g-4 g^Q]-^ gmall fields are well below the Paschen-Back regime, which is reached 
around 10"'^ - 10"^ T for the npi/2 - 7^^3/2 and 72^3/2 - nd^/2 fine structure splittings at 
= 50 [35^. Hence, both (. and j are good quantum numbers, and additionally, any 
splitting of the different mj levels of a given j state are much smaller than the pulsed 
laser linewidth. In other words, we can be certain which n^j state we are exciting, and 
there is no ^-mixing of states with different lifetimes. 

The second issue is that there is a small residual electric field in the interaction 
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region, which could lead to Stark mixing of different i states. The electric field is due 
primarily to the bias voltage on the MCP, which (at the time of these experiments) 
was not housed in a grounded enclosure, though it was separated from the interaction 
region by the meshes described in section I3.2[ We have used this same apparatus to 
study ultra-cold plasmas (similar to those described in 2l|] and 50|). In the plasma 
experiments we found it necessary to apply a nulling field of no more than 1 V/cm 
to maximize the plasma lifetime. However, this field was not corrected in the lifetime 
measurements described in this paper, and we therefore estimate that the Rydberg 
atoms experienced a static electric field of < 1 V/cm. The Stark interaction mixes 
states with orbital angular momentum i with states having i' = i±l (which themselves 
mix with states ± 1). However, (or, in atoms with large fine structure splittings 
like Rb, rrij) remains a good quantum number. This mixing obviously has the potential 
to impact our lifetime measurements, given the generally much longer lifetimes of the 
np states than of the ns and nd states (whose lifetimes are similar) at comparable UqQ. 

To estimate the extent to which our measurements are an average over Stark-mixed 
ni states, we have calculated the amount of mixing that could result from an electric 
field using a program that diagonalizes the Stark energy matrix 5l|] . For a given electric 
field, the mixing is stronger as n increases, so we calculated the mixing at the highest 
rZgff for any state whose lifetime was measured in our experiments (44^5/2, corresponding 
to rigfj- ~ 43). Additionally, we assumed an electric field of 2 V/cm, i.e., twice the field 
that was present in our experiments. We found that in this regime, the admixture of 
states with a different orbital angular momentum quantum number i was negligible. For 
instance, the 44^5/2 \mj\ = 1/2 state is 1.57% (probability) i = I, 97.10% i = 2, and 
1.29% i = 3, while the 46si/2 \mj\ = 1/2 state is 98.95% £ = 0, and 1.05% i = 1. We 
looked at all the possible \mj\ eigenstates in the 46si/2, 45pi/2, 45p3/2, 44^3/2 and 44^5/2 
states in a 2 V/cm field, and found that none had less than 97.1% (by probability) of its 
zero- field i eigenstate content, and the maximum amount of a state with a different i 
content was the 1.57% p-character in the 44^5/2 \mj\ = 1/2 state. Such small admixtures 
do not impact the measured state lifetimes at a level that can be distinguished in this 
experiment. 

We therefore conclude that the residual electric and magnetic fields have no 
impact on our lifetime measurements. (The < 5 V/cm voltage pulse which is applied 
simultaneously with the laser when np states are to be excited decays to zero over ^ 1 /xs. 
Hence, the pulse switches off slowly enough that the Stark state excited adiabatically 
couples to £ = 1 in zero field and introduces no ^-mixing [52].) 

Our measurements roughly correspond to those made using_^ a slightly different 



method by a group at Universidade de Sao Paulo, Brazil [27|; |28j. Their full set of 
results for ns (28 < n < 44) and nd (27 < n < 42) states is given in 28|] and 29||, and 
their results for the np (33 < n < 42) in |22] (their laser excitation method prevented 



were 



isolation of the different j-states for np and nd). The np lifetimes given in [27 
erroneously reported to be twice the actual values, as they pointed out in the reference 
section of their subsequent paper [2^, and their np lifetime results were afterwards 
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Figure 5. Our results for the radiative lifetimes, of (a), the Rb nsi/2 

(28 < n < 45) states; (b), the np3/2 (34 < n < 44) states; and (c), the 
(29 < n < 44) states, at a nominal radiation temperature of T = 300 K, plotted 
versus the effective principal quantum number, n^Q- The graphs also show theoretical 
predictions based on ([T]), ([5]), ([2]), and ([3]), using parameter values reported in a number 
of different papers. All the models use the BBR rates of Beterov et al (equation (14) 

in 0); the solid red line ( ) uses values for rs and e from Beterov et aZin ([2]); the 

purple dashed ( ) line uses Gounand's values for ts and £^.^1; the blue dash- 
dot-dash ( — • — ) Hue uses the values of He et al for ao....a3 in ([3]) 5]. The predictions 
using the values of He et al for ts and e in ([2]) are essentially indistinguishable from 
those found using their values for ao....a3 in ([3]), except for the curve, which is 

shown as a green dotted ( ) line. 
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Table 2. Our experimentally determined radiative lifetimes for the nsi/2, and 
Rydberg states of Rb at a temperature of 300 K. For the nd state lifetimes 
denoted by the asterisk, the separation of the nd j = 3/2 and j = 5/2 states is less 
than the bandwidth of the pulsed laser used to populate the target state at i = 0. The 
measured lifetime is therefore likely subject to some degree of averaging over the two 
j states. 
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39 


34.3 ± 1 


5 
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40 
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52.6 ± 1.3 
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60.2 ±2.8 


38 


9 ±3.0* 


44 


47.0 ± 1 


7 


64.2 ±2.6 


42 


0±3.2* 


45 


50.8 ±2 


4 









reported correctly in 2^. (This group's experimental results for the ns, np, and nd 
lifetimes are also reproduced in full in figure 6 of the paper by Beterov et al [6|].) In 
general, our s- and d-st&te results are higher than the Sao Paulo group's measurements 
by an average of 7-8%. While some values differ by 20-30% for low n, there are negative 
as well as positive differences between the data sets. Furthermore, our p-state results 
are all lower than the previous measurements by between 5% and 16%. However, the 
difference between the two data sets is not large, and almost all the lifetimes agree within 
three times the combined experimental uncertainty. We conclude that while there seem 
to be systematic differences between our results and those of the Sao Paulo group, 
these differences are comparable with the measurement uncertainties. Tretyakov et al 
have also measured the lifetime of the 37p state of Rb using a MOT source, obtaining 



r = 38 ± 3 /is, which is lower than our value of 43.5 ± 1.9 /is 31 



In a recent paper reporting the first experimental observation of a new class of 
Rydberg molecules, lifetime measurements were used to distinguish signals from atomic 
and molecular species. Specifically, Bendkowsky et al observed ^^Rb2 dimers comprising 
one atom in the 5si/2 ground state, and the other atom in a nsi/2 Rydberg state, 
where 34 < n < 40 [2(J\. The molecular resonances were observed as satellites to 



Radiative lifetime measurements of rubidium Rydberg states 



19 



atomic resonances when a tunable 480 nm cw laser was scanned over the 5j93/2 — ^ ^■§1/2 
transition, and the resonances were observed using SFI followed by detection of the 
resulting ions with an MCP. The atomic sample had a temperature of 3.5 /iK and a 
density of 1.5 x 10^^ cm~^, i.e., similar conditions to those used to create a quantum- 
degenerate gas of Rb. Bendkowsky et al report a number of experimental observations 
that strongly justify their claim that they observe new molecular species, including 
measurements of lifetimes for the different spectral features. While their method was 
able to distinguish different decay times for the atomic and molecular states, their 
measured lifetimes of the atomic states are very different from our own (and those of 
the other experimental works discussed in the previous paragraph), ranging from 65 
± 9 /is for 35s to 57 ± 5 /xs for 37s. These differences are surprising, even when 



it is considered that in 20|, the detection of positive ions meant that the observed 



signal was insensitive to redistribution by BBR so that many Ryd berg states likely 



contributed to it. On the other hand, the experiments reported in [20| were carried 
out with MOT atom densities ~ 10^ times larger than in our own study and those 



26 



27 



28 



29 



30 



31 



32l |. Indeed, it is unclear that an atomic property 



reported in 

like a radiative lifetime can be measured under these conditions, give n the occurrence 
of many-body effects at much lower densities [l3|; Q; M, H; M ; [ISJ. It is, however, 
curious that in this instance, any such effects (or other mechanism) cause longer state 
lifetimes to be observed for dense samples than for individual atoms. 



4-2. Comparison with theory 

We considered three different means of comparison between our data and existing 
theoretical results. First, we estimated which theoretical result was most compatible 
with our data set. Second, we combined our data with data from previous \ow-n 
measurements to find experimental values for Ts and e in ([2]). Finally, we analyzed 
the sensitivity of these parameters to our imperfect knowledge of the BBR temperature. 



4-2.1. Least-squares comparison with theory We synthesized theoretical lifetime values 
based on ([1]), ([5]), ([2]), and ([3]), using parameter values reported in (we used 

the values denoted by the asterisk (*) in table [H for the np3/2 lifetime predictions of 
Beterov et al [6]). All the models use ([5]) for the BBR rates from j^. These synthesized 
lifetime values are shown as smooth curves in figure [5l The behaviour of the zero-Kelvin 
lifetimes in the work of Gounand and of Beterov et al jsl; @] is of the form of (El), while 
He et al also use this form, as well as that of ([3]) [5|. (The predictions using values for 
Ts and e from He et al in ([2]) are essentially indistinguishable from those found using 
their values for ao-.-.a^ in ([3]), except for the nd5/2 states.) 

We used a variant of the "chi-squared" statistical parameter to compare our results 
with the theoretical curves. Specifically, we calculated the values for 

N 

,2 _ (Exp Jrh\2 I ( Exp\2 /n\ 

i=l 
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where t^^^ are our measured values, the theoretical values, and a^^^ our 

experimental uncertainties. These parameters were calculated for each of the nsi/2, 
^P3/2; and series. For each angular momentum series, is the number of different 
state lifetimes that were measured. 

For the nsi/2 and np3/2 states, is consistently lower for the results of Beterov et 
al , although the predictions of He et al are almost as good, the values predicted 
by ([3]) being slightly better than those using ([2]) (We have been advised by one of 
the authors of [6|] that the and e for the np3/2 states are in error 36|. The corrected 
values are given in table [1] of the present paper, and are to appear in an erratum to 
[gI. Using the revised values, the lifetime predictions for the states almost exactly 
match those of He et al using either of ([2]) or ([3]).) For the nd^i2 states, the results of He 
et al using ([3]) agree with the data slightly better than those of Beterov et al , with the 
prediction of He et al from ([2]) being worse than either. For all of the measurements, 
Gounand's predictions jl] have a much poorer agreement than either of He et al results 
or that of Beterov et al . 

We conclude that our data set is unable to definitively distinguish which of the 
works of Beterov et al or He et al gives better parameters Tg and e in ([2]), nor whether 
([2]) or ([3]) given by He et al is more accurate. However, it seems that the accuracy of 
the theoretical calculations of both of these theoretical papers is better than that of 
Gounand jl]. 

4..2.2. Experimental dependence of tq on n^jj We initially tried fitting our data to ([2]), 
assuming a radiation temperature of 300 K and the equation for the BBR depopulation 
rate, ([5]), from Beterov et al . That is, we equated our experimental lifetimes, r^^(T), 
with l/Fgff (300 K) in ([1]) assuming that ([5]) was correct, to find Fq and hence 
experimental values for tq (designated Tq^^): 

We then used to extract experimental values for Tg and e (assuming that T = 300 
K), using a least-squares fit. However, since our measurements were taken at relatively 
high n, there is a large uncertainty in the values obtained, ~ ±100% for Ts and ~ ±10% 
for e. We concluded that using only our own lifetime results was insufficient to make a 
meaningful comparison with theory. 

To obtain more accurate values for Tq and e, we added data points from previous 
measurements of Rydberg state lifetimes in Kb. Specifically, we included values for the 
following states: nsi/2 with n = 10, 11, 12, 14, 16, 17, 18; np3/2 with n = 10, 11, 12, 14, 
17, 22; and nd^/2 with n = 10, 11, 12, 13, 15, 18. These experimental values are the most 
precise measuremen ts g iven in table V of Theodosiou's paper 0], and were originally 



reported in |38|; |4l|; |42|]. Using these data points, as well as our own, we plotted the 



F^"^^ values obtained using ([8]) versus UqQ, and these data are presented in figure [61 
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Table 3. Values of the parameters rg and e in ([2]) obtained by fitting this equation to 
the data given in figure [6l 



State 




e 




(ns) 




nsi/2 


f.4±0.1 


2.99 ±0.03 




3.5 ±0.4 


2.90 ±0.03 




1.8 ±0.3 


2.84 ±0.04 



We excluded previous data for n < 10 because the accuracy with which 



reproduces the numerical result for Tggj^ is unspecified below 



n = 10 (6|. When 

obtaining the zero-Kelvin decay rates shown in figure [6], we assumed a BBR spectral 
temperature of exactly 300 K for our own experiment. In fact, this may not be the case, 
an issue which is considered below in section 14.2.31 With regard to the previous results 
shown in figure El we have used the estimated BBR temperatures given in those papers. 
Specifically, in |42l |. the Rb was contained in a 320 K cell; in 38|, the sample was at 460 
K; and in 



41 



the Rb was maintained at 520 K. We have assumed these temperatures 
to be exact, i.e., the uncertainty in the F^^^ values is entirely due to the experimental 
uncertainty of the lifetime measurements, with no contribution from rBgj^(T). 

We fitted the data presented in figure [6] to ([2]) to obtain experimental values for 
Ts and e for the nsi/2, np3/2, and 71^5/2 states of Rb (see table [3]). Comparing these 
experimental results to the theoretical values given in table JT], it appears that the 
dependence of Tq^^ on n^Q is best reproduced by He et al |5|, rather than by that 
of Beterov et al (in which the parameter e shows essentially no dependence on t). The 
three theoretical papers substantially agree on the Tg values for the nsi/2 and 

npj series, and our value of Tq for nsi/2 is consistent with this finding, although our 
value for npj is higher than the theoretical values by ~30%. For the ndj data, the rg 
value best agrees (within 10%) with the result from He et al , being 15% lower than 
Gounand's and 70% higher than that of Beterov et al . 

It is unclear whether any differences of the experimental Tg values from those given 
in any of the three theoretical papers is meaningful, given the 10% uncertainty of the 
experimental values. On the other hand, a preliminary consideration of the experimental 
e values obtained from the extended data set seems to favor the predictions of He et al 
, rather than those of Gounand or Beterov et al . However, we believe that any such 
comparison must be treated with caution, given the uncertainty in our knowledge of the 
BBR temperature, which we discuss below. 

Nascimento et al H found rg = 1.43 ±0.05 ns and e = 2.94 ±0.03 for the Rb 
states in the range 28 < n < 44, and rg = 1.90 ± 0.07 ns and e = 2.83 ± 0.04 for the 



Rb ndj states in the range 27 < n < 42 (these values are also reported in [29 




results agree quite well with our findings shown in table [3l and with theory 
However, as one of us has observed 33|], any agreement of the results in 28| with ours, 
or with the theoretical values, is apparently accidental, since Nascimento et al used the 
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Figure 6. (a) Log-log plot of an extended set of rubidium Rydberg state zero-Kelvin 
decay rate measurements, T^^^, comprising our results and those of 
nsi/2 (full red diamonds, ♦), np3/2 (open green circles, O and nd^^2 (blue crosses, 
+) states with n > 10. The range of ti^q values is 6 < ti^q < 50, and the range of T^^^ 
values is 4 X 10'^ s^^ < Tq^^ < 4 x 10^ s~^. The previous results are the most precise 
values given in table V of Q that satisfy the criterion that n > 10. The zero-Kelvin 
decay rates are obtained from the data using ([S]), assuming that the calculations of 
Beterov and co-workers for the BBR depopulation rates are correct. The straight lines 
are fits of the data to i.e., we fit the data to the equation Tq^^ = (l/''s) "•gjj- The 

full line ( ) is the fit to the ns 1/2 data, the dotted line ( ) to the np^^2 data, 

and the dashed line ( ) to the 71.^5/2 data. The dotted box shows our data, and is 

expanded in graph (b) below, (b) Inset of data presented in (a) between 24 < rigg- < 44 
showing only our results as a linear-linear plot. The y-axis is calibrated in units of 
1000 s~^. This plot also shows the fit lines to the entire data set shown in (a). 
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approximate equation (jl]) for BBR depopulation rates to find the zero-Kelvin lifetimes 
from their experimental data. As noted above, the approximate equation overestimates 
the BBR depopulation rates by between 20% and 40% at n = 30 in Kb. (A very recent 
experimental paper on cesium Rydberg state lifetimes also uses the approximation for 



the BBR rates when comparing results with theory 321], as does 291] when comparing 
the Ts and e values for Rb np states with theory.) 

4-2.3. Effect of BBR temperature The ability of BBR to significantly affect state 
lifetimes was investigated in several studies in the late 1970's M- We have 

not attempted to measure or control the BBR temperature with any high degree of 
precision in our experiment, and this approach is similar to other recent Rydberg lifetime 



experiments 26|; |27|; l28|; l29|; l30|; l3l|; 1321]. Our lab is stabilized at 295 K; however. 



numerous pieces of equipment become warm during use, and, as pointed out by Spencer 



et al [id], it is very difficult to shield any experiment from BBR from external sources. 
Most of the BBR absorbed by atoms in our apparatus will come from the walls of 
the vacuum chamber, and will have a radiation temperature close to 300 K. However, 
the atoms will also be illuminated by BBR from the Rb side-arm, and this source was 
heated to between 330 and 350 K. In other words, there is no single temperature that 
characterizes the BBR field, since the vacuum chamber is not in thermal equilibrium. 

To investigate the potential impact on the values of Fq and tq that we inferred from 
our experimental data, we performed a simple analysis using expressions from Beterov 
et al for tq ((Ej) and the values for Ts and e given at the bottom of table [I]) and Tggj^ 
(([5]) and the values for A, B, C, and D given in table I in j^). We make the gross 
assumption that ([2]) and ([5]), when used with parameter values from Beterov et al , 
exactly predict the result of an experiment carried out on atoms illuminated with BBR 
with characteristic temperature T. In this situation, the lifetime values measured, t^^, 
would be 

-E^, = ^+^lun (9) 

On the other hand, if these values of t^^^ are used to infer values for the zero-Kelvin 
lifetime, Tq"^, assuming that the BBR temperature is 300 K, we have 
1 1 



E.,^^,- W + WSOOK). (10) 

Obviously, if T = 300 K, fitting the values of Tq"-^ to ([2]) would yield the original values 
for Tg and e. The question we want to investigate is, what happens to the values of 
these parameters obtained if T 7^ 300 K? We consider two situations: T = 270 K, and 
T = 330K, and made another gross assumption, namely that the "experimental" values 
of Tq"-^ have uncertainties that are a fixed percentage of their absolute values. (If the 
data points are unweighted, the low-n lifetimes have a much smaller impact on the fit 
than the high-n lifetimes.) Finally, we assume that we "measured" lifetimes of only the 
states shown in figure El 
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We did this analysis for the ns, np, and nd series, but we will discuss only the 
np results. Because the np states have longer zero-Kelvin lifetimes, BBR has a bigger 
impact on the effective np lifetimes, and therefore a wrong temperature assignment has 
a bigger impact on the inferred zero-Kelvin lifetime. We found that while a difference 
of 10% between the actual value of T and its assumed value (300 K) changed e by only 
3%, from approximately 3.0 to 2.9 or 3.1, the value of Ts changed by about 15%, from 
2.5 ns to 2.1 ns or 2.9 ns. (The range of these values does not correspond to the Tq value 
for the np3/2 states given in table [3] since we are starting with theoretical values from 
[gI, not our own experimental data.) As expected, the effect on the ns and nd results 
was smaller. 

We conclude that the ability of experimental measurements of Rydberg state 
lifetimes to inform the results of theoretical calculations is substantially dependent on 
accurate knowledge of the BBR field temperature. This is an issue that can only be 
resolved if BBR depopulation rates of the Rydberg states can be measured under the 
same experimental conditions as those of the lifetime measurements. It is unclear, 
however, whether the BBR transition rates can be measured with sufficient precision to 
yield zero-Kelvin lifetimes with small enough uncertainties for useful comparison with 
theory. The two studies that have specifically addressed this question, those of Spencer 



et al lOj and Galvez et al 12|, did not measure the BBR transfer rates explicitly. 



Rather, they showed that their models (based on theoretical rates) were consistent with 



the observed time-evolution of population over a range of Rydberg states 12| , or that the 



observed variation of the radiative lifetime of a state with temperature was in agreement 



with calculations 10 



5. Conclusion 

We have measured the radiative lifetimes of high-?7, Rydberg states of rubidium using 
cold atoms in a MOT. The slow velocities of these atoms allow them to be observed 
for long enough that lifetimes of states up to at least n = 50 may be measured. Our 
experiment has a significant improvement over earlier work that reduces the impact of a 
possible source of systematic error. Specifically, we use mm-waves to "tag" a particular 
state so that the detected signal is an accurate mirror of the target state population. 
This eliminates ambiguities in the signal caused by BBR stimulated population transfer 
to nearby Rydberg states. 



Our results are consistent with other recent experimental studies 27|; l28|; 1291], and 



; experim( 

also agree well with several theoretical paper. 113. However, we hnd that the 
ability of our data to distinguish the merits of the three different theoretical studies 
is limited by a lack of experimental data of BBR depopulation rates in rubidium, and 
by our imperfect knowledge of the BBR spectral temperature. On the other hand, we 
have used available theoretical information on BBR rates to obtain the highest possible 
accuracy for our values for the zero-Kelvin lifetimes of Rb Rydberg states. 



REFERENCES 



25 



Acknowledgments 

We acknowledge many conversations with Mike Noel, and especially Tom Gallagher, 
Wenhui Li, and Paul Tanner, who gave invaluable assistance with the design of the mm- 
wave spectrometer and the narrow-bandwidth pulsed laser system. Charlie Conover 
provided much assistance on calculating Stark-mixing effects due to residual electric 
fields, as well as in many other matters. We appreciate the willingness of Phillip Gould 
to share with us plans of the mount for the high-transparency field ionization meshes. 
Ryan Jennerich and Greg Foltz assisted with the construction of the MOT used in this 
experiment. This work has been supported by NSF (grant numbers OlA-9873763, PHY- 
0140430, and PHY-0652842) and by Colby College via the Division of Natural Sciences 
Grant Program. Tamas Juhasz and Cristian Vesa acknowledge support from Colby 
College in the form of summer research assistantships, and Roy Wilson is grateful to the 
NSF for a supplemental summer research assistantship (grant number PHYS-0734190). 

References 

[1] Hansch T W 1972 Appl. Opt. 11 895 

[2] Littman M G and Metcalf H J 1978 Appl. Opt. 17 2224 

[3] Gounand F 1979 J. Physique (Pans) 40 457 

[4] Theodosiou C E 1984 Phys. Rev. A 30 2881 

[5] He X, Li B, Chen A and Zhang C 1990 J. Phys. B: At. Mol. Phys. 23 661 

[6] Beterov I I, Ryabtsev I I, Tretyakov D B and Entin V M 2009 Phys. Rev. A 79 
052504 

[7] Gallagher T F and Cooke W E 1979 Phys. Rev. Lett. 42 835 

[8] Ducas T W, Spencer W P, Vaidyanathan A G, Hamilton W H and Kleppner D 
1979 Appl. Phys. Lett. 35 382 

[9] Belting E J, Hildebrandt G F, Kellert F G, Foltz G W, Smith K A, Dunning F B 
and Stebbings R F 1979 J. Chem. Phys. 70 3551 

[10] Spencer W P, Vaidyanathan A G, Kleppner D and Ducas T W 1982 Phys. Rev. A 
25 380 

[11] Farley J W and Wing W H 1981 Phys. Rev. A 23 2397 

[12] Galvez E J, Lewis J R, Chaudhuri B, Rasweiler J J, Latvakoski H, De Zela F, 
Massoni E and Castillo H 1995 Phys. Rev. A 51 4010 

[13] Anderson W R, Veale J R and Gallagher T F 1998 Phys. Rev. Lett. 80 249 

[14] Mourachko 1, Comparat D, dc Tomasi F, Fioretti A, Nosbaum P, Akulin V M and 
Fillet P 1998 Phys. Rev. Lett. 80 253 

[15] Tong D, Farooqi S M, Stanojevic J, Krishnan S, Zhang Y P, Cote R, Eyler E E 
and Gould P L 2004 Phys. Rev. Lett. 93 063001 



REFERENCES 26 

[16] Singer K. Rcctz-Lamour M, Amthor T, Marcassa L G and Weidemiiller M 2004 
Phys. Rev. Lett. 93 163001 

[17] Cubel-Liebisch T, Reinhard A, Berman P R and Raithel G 2005 Phys. Rev. Lett. 
95 253002 

[18] Vogt T, Viteau M, Zhao J, Chotia A, Comparat D and Fillet P 2006 Phys. Rev. 
Lett. 97 083003 

[19] Stanojevic J, Cote R, Tong D, Farooqi S M, Eyler E E and Gould P L 2006 Eur. 
J. Phys. 40 3 

[20] Bendkowsky V, Butscher B, Nipper J, Shaffer J P, Low R and Pfau T 2009 Nature 
458 1005 

[21] Li W, Noel M W, Robinson M P, Tanner P J, Gallagher T F, Comparat D, 
Laburthe-Toha B, Vogt T, Zahzam N, Vanhaecke N, Fillet F and Tate D A 2004 
Phys. Rev. A 70 042713 
[22] Li W, Tanner F J and Gallagher T F 2005 Phys. Rev. Lett. 94 173001 
[23] Li W, Tanner F J, Jamil Y and Gallagher T F 2006 Eur. J. Phys. D 40 27 
[24] Li W, Mourachko 1, Noel M W and Gallagher T F 2003 Phys. Rev. A 67 052502 

[25] Han J, Jamil Y, Norum D V L, Tanner F J and Gallagher T F 2006 Phys. Rev. A 
74 054502 

[26] Magalhaes K M F, de Oliveira A L, Zanon R A D S, Bagnato V S and Marcassa 
L G 2000 Opt. Commun. 184 385 

[27] de Oliveira A L, Mancini M W, Bagnato V S and Marcassa L G 2002 Phys. Rev. 
A 65 031401 (R) 

[28] Nascimento V A, Caliri L L, de Oliveira A L, Bagnato V S and Marcassa L G 2006 
Phys. Rev. A 74 054501 

[29] Marcassa L G 2009 Phys. Scr. T134 014011 

[30] Gabbanini C 2006 Spectrochim. Acta, Part B 61 196 

[31] Tretyakov D B, Beterov I I, Entin V M, Ryabtsev I I and Chapovsky P L 2009 
JETP 108 374 

[32] Feng Z G, Zhang L J, Zhao J M, Li C Y and Jia S T 2009 J. Phys. B: At. Mol. 
Phys. 42 145303 

[33] Tate D A 2007 Phys. Rev. A 75 066502 

[34] Cahri L L and Marcassa L G 2007 Phys. Rev. A 75 066503 

[35] Gallagher T F 1994 Rydberg Atoms (Cambridge, UK: Cambridge University Press) 

[36] Beterov I I 2009 Private communication 

[37] Cooke W E and Gallagher T F 1980 Phys. Rev. A 21 588 

[38] Gounand F, Fournier P, Cuvellier J and Berlande J 1976 Phys. Lett. A 59 23 

[39] Lundberg H and Svanberg S 1976 Phys. Lett. A 56 31 

[40] Hugon M, Gounand F and Fournier P R 1978 J. Phys. B: At. Mol. Phys. 11 L605 



REFERENCES 



27 



[41] Gounand F, Hugon M and Fournier P 1980 J. Physique (Paris) 41 119 

[42] Marek J and Munster P 1980 J. Phys. B: At. Mol. Phys. 13 1731 

[43] Jeys T H, Foltz G W, Smith K A, Beiting E J, Kellert F G, Dunning F B and 
Stebbings R F 1980 Phys. Rev. Lett. 44 390 

[44] Monroe C, Swann W, Robinson H and Wieman C 1990 Phys. Rev. Lett. 65 1571 

[45] Arnold A S, Wilson J S and Boshier M G 1998 Rev. Sci. Instrum. 69 1236 

[46] Arimondo E, Inguscio M and Violino P 1977 Rev. Mod. Phys. 49 31 

[47] Lee S A, Helmcke J, Hall J L and Stoicheff B P 1978 Opt. Lett. 3 141 

[48] Ye J, Swartz S, Jungner P and Hall J L 1996 Opt. Lett. 21 1280 

[49] Thibault C, Touchard F, Biittgenbach S, Klapisch R, de Saint Simon M, Duong 
H T, Jacquinot P, Juncar P, Liberman S, Pillet P, Pinard J, Vialle J L, Pesnelle A 
and Huber G 1981 Phys. Rev. C 23 2720 

[50] Killian T C, Kulin S, Bergeson S D, Orozco L A, Orzel C and Rolston S L 1999 
Phys. Rev. Lett. 83 4776 

[51] Zimmerman M L, Littman M G, Kash M M and Kleppner D 1979 Phys. Rev. A 
20 2251 

[52] Freeman R R and Kleppner D 1976 Phys. Rev. A 14 1614 



